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We describe an entanglement witness for iV-qubit mixed states based on the properties of Ap- 
point correlation functions. Depending on the degree of violation, this witness can guarantee that 
no more than M qubits are separable from the rest of the state for any M < N, or that there is some 
genuine M-party or greater multipartite entanglement present. We illustrate the use our criterion 
by investigating the existence of entanglement in thermal stabilizer states, where we demonstrate 
that the witness is capable of witnessing bound-entangled states. Intriguingly, this entanglement 
can be shown to persist in the thermodynamic limit at arbitrary temperature. 



I. INTRODUCTION 

Recently, there has been growing interest in the en- 
tanglement of large many-body systems. The reasons for 
this are many and varied. From the quantum information 
perspective, quantum computation intrinsically involves 
the use of large arrays of entangled qubits, including the 
one-way model of computation [T] where the initial, clus- 
ter state, of the system is considered a multipartite en- 
tangled resource for computation. This model for compu- 
tation has the potential to enable quantum information 
processing in many physical systems which were other- 
wise considered infeasible. Generating the cluster state 
and finding signatures with which to recognise its suc- 
cessful generation along with protecting and recovering 
it from the effects of noise are thus matters of great im- 
portance. Of more immediate relevance to condensed 
matter theory and experiments is the study of, for exam- 
ple, how thermodynamic parameters of many-body sys- 
tems as well as phase transitions can be directly related 
to entanglement [3JH]. However, it is still not clear if this 
relation always holds and if it is of any real significance 
to our understanding of the physics of large quantum 
systems. 

To be able to better investigate such properties both 
theoretically and experimentally, it is interesting to de- 
velop techniques for entanglement detection in many- 
qubit systems. The most straightforward approach is 
via the so-called entanglement witnesses. Entanglement 
witnesses are Hcrmitian operators such that their mean 
value with an arbitrary separable state is always less than 
one. Thus, if the observed value exceeds one, we are cer- 
tain that a given state is entangled. The significance 
of the entanglement witness approach is that, in princi- 
ple, it can be experimentally implemented and that one 
can find entanglement witnesses tailored to detect multi- 
partite entanglement [H |7] as well as entanglement 
strong enough to violate some Bell inequalities. It was 
shown in Ref. [5] that one can perform certain computa- 
tional tasks more efficiently than using classical resources 
only if entanglement violating Bell inequalities is avail- 
able as a resource. 



There is little literature dealing with entanglement wit- 
nesses for many qubit systems capable of detecting gen- 
uine multipartite entanglement [10J [TTJ [12] . The usual 
approach (see, for instance, Ref. [7]) is to find an en- 
tanglement witness tailored for a specific state and then 
show that it can be measured locally. Usually such wit- 
nesses are capable of detecting entanglement only in the 
neighbourhood of the state for which they have been op- 
timized. In this paper we derive a witness that is capable 
of detecting multipartite entanglement for a wide set of 
multi qubit states. Moreover, the number of different set- 
tings required to measure our witness does not have to in- 
crease exponentially with the number of qubits (although 
generally it does), which is often the case [THl HU fT2"] . 
with the notable exceptions of Refs. [3 [8]. 

In this paper, we find a family of entanglement wit- 
nesses for N qubits, the mean value of which solely de- 
pends on iV-point correlation functions. We derive a rel- 
atively simple lower bound for the witness and demon- 
strate its calculation with the examples of thermal sta- 
bilizer states and the thermal single-excitation Bose- 
Hubbard model to obtain a range of temperatures for 
which entanglement exists. We also relate the lower 
bound to the violation of Bell inequalities from the so- 
called Werner- Wolf, Zukowski-Bruckner (WWZB) fam- 
ily [T5]. We show that depending on the different level 
(strata) of violation of the inequalities, we can assign 
lower bounds to the degree of entanglement present. 



II. WITNESSING SEPARABILITY AND 
MULTIPARTITE ENTANGLEMENT 

In |14j , a convenient parametrization of all two-setting 
Bell inequalities of the WWZB type was developed by 
considering the Hcrmitian operator 

fce{o,i} N 

subject to a variety of constraints on the coefficients b^. 
The operators = \Gi)(Gi\ are orthogonal projectors 



2 



on the generalized GHZ states |G0 = i(|jfe) iof-^fc)). 

For convenience, we can express af N \k) = \k'). We will 
now investigate the properties of W devoid from the re- 
strictions imposed in [14] , except that we shall require W 
to be an entanglement witness. 



A. Conditions for Full Separability 

The family of operators UWU\ where U = ]J^ =1 U {n) 
and CTtO is an arbitrary SU (2) transformation on qubit 
I, becomes a family of entanglement witnesses if 



(2) 



B. Entanglement Witnesses 



In order to detect entanglement, we have to calculate 
Tr(TFp) =m&xJ2 b k Re ($\UpU t \k')) 

U k 

subject to the constraint ([2]). We shall now consider two 
specific cases of entanglement witness W, by selecting 
two specific sets of b^. The first, Wa, is the strongest 
witness of this class, and is selected by finding the value 
of fco which maximizes max^ |Re((/c|C/p[/^|fc'))|, setting 
br to 2^, and all others to 0. 

fco 

In general, it is likely that the maximization over k 
will be difficult to do, so we choose another witness, Wb, 
which will give a lower bound to this value. We select 



since, if this condition holds, one has 



(ip sep \uw&\*p s 



< 1 



sep " " I Y sep/ 

for an arbitrary pure and fully separable state 

N 



(3) 



\^sep)=6dV^\d). 



In the above formula, is an SU(2) transformation 
which can be parametrized by = e - ien j- a ; where n 
is a unit vector. To prove Eqn. (jjjl, firstly observe that 
the unitarics can be absorbed in so that we 

only have to show that 



max(^ sep |H / |V' se p) 

Wsep) 



max 

k 



^2 b k Re ( (^ep | k) (k'\tp sep ) 



< 1, 



(4) 



where we've made use of the expansion W = 
\ Eg h $\ + W) (k\) ■ The properties of im- 
pose that 

Re((^ep|0)(0> sep )) = Rem(0|FWt|Q){l|V-W|Q)] , 



where we have also absorbed the a x rotations from k to 
into the Therefore, the maximum over n and 9 

of Eqn. Q yields 



max{ip aep \W\ip sep ) 



2 N 



(5) 



with values n z = n x = I and = 0. This is not greater 
than 1 when Eqn. ^ holds. Proving this for pure states 
is sufficient since the convexity of the mixed separable 
states implies that the optimum will be given by a pure 
state. 



where 



H = J2(-lf T co S (*\t\)T r . 



(6) 



The number T f = Tt{Uap^ p) is the average value of spin 
measurements along the directions given by O n x, O n y 
[n = 1, . . . , N), where O n is an orthogonal representa- 
tion of U n . We see that the only relevant iV-point corre- 
lation functions are those for which |Z|, i.e., the number 
of a y in a-j- before the local unitary operation U, is an 
even number. The choice of b^ is clearly sub-optimal, 
so Ti(Wap) > Tt(Wbp)- The reason for this particular 
choice is that 



Tr(Wsp) = max 
u 



E 

\l£even 



(7) 



as will be proved in the Appendix. 

Our starting point was motivated by the fact that it 
was shown in Ref . [2] that the family of operators IJWW 
subject to the constraint J^k b % ~ ^ N an( ^ some com- 
plicated additional constraints on the signs of the &g's 
coincides with the family of all two-setting Bell inequal- 
ities of the WWZB type [15] . Therefore, if for some 
state p one maximizes Tr (W^p) only under the constraint 
Efc^| = 2 N , one gets the upper bound for violation 
of two-setting WWZB inequalities. If this maximum is 
larger than one, nothing conclusive can be said about vi- 
olation of two-setting WWZB inequalities for the state p. 
However, if Tr(W^p) < 1, one concludes that the state 
p cannot violate any of the two-setting inequalities from 
the WWZB family. 

Specifically, selecting — 2 N A fc / \/J2k ^fc yields the 
maximum under the constraint Efc^l = ^ ' an d gives 



W c 



max 

u ' 



iGeven 



rp2 



(8) 
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Thus, if Wc < 1, which is equivalent to the condition 
Wb < 1, one cannot violate two-setting Bell inequalities 
from the WWZB family. This is in agreement with the 
necessary condition given in |13j for not violating two- 
setting WWZB inequalities. Interestingly Wb > 1 is 
a sufficient condition for violation of the multi-setting 
WWZB Bell inequalities presented in Ref. [TS] . 

C. Wa as a witness of partial separability 

Let us examine the witness Wa more carefully detail- 
ing its effect on pure states p = \i/)){ip\. 

Tr{W A p) = 2 A 'maxRe((fc|J7|V>)(^ t |fc')) 

k,U 

= 2 Ar maxRe ((0\ 8N U\il>)(il>\ff\l)* N ) . 

We can now split \ip) into a part \tp) acting on M qubits, 
and \ipsep): a m lly separable state on the other N — M 
qubits. The witness does not change under permutations, 
so without loss of generality, we can take the qubits of 
|^) to be the first M qubits. Having already established 
that 

m^M(0\ 9N - M U\^){^\t^\l) 9N - M ) = 

we see that the maximum value of Tr (WaP) where 
no more than M qubits are entangled is given by 
Tr(W^4|V')(V , l)> our entanglement witness on M qubits. 
Hence, the largest possible value is 2 M_1 , where — 
(|0)® M + |1)® M )/V2. We conclude that if Tr(W A p) > 
2 M ~ 1 , fewer than N — M qubits are separable. Again, 
by convexity, the results also apply to mixed states, even 
though we only performed the calculation for pure states. 
So, our entanglement witness witnesses not only the fact 
that there is some entanglement, but witnesses that at 
least M + 1 qubits are entangled, although this cur- 
rently makes no statement about the type of entangle- 
ment present. 

D. Wa as a witness of multipartite entanglement 

Given that we can witness the fact that many qubits 
are entangled, it would also be interesting if we can wit- 
ness different types of multipartite entanglement. The 
first step in this process is to determine if M-party (or 
greater) multipartite entanglement is involved. Again, 
we consider W A acting on pure states. Since we know 
that GHZ states maximize the value of the witness (by 
design), then if M -partite entanglement is involved (as- 
suming N/M is an integer), the maximum violation of 
Wa must be given by = \/p)® N / AI , where \tp) — 
(|0)® M + \1)® M )/V2. We find that 

Ty(W aP ) = Tr(W A $)$\) N/M = 2 N - N ' M . 



Thus, if Tr(W A p) > 2 N ~ N / M , the multipartite entangle- 
ment that is involved must be at least (M + l)-partite. 
If N/M is not an integer, the maximal violation is given 
by \N/M\ copies of an M-qubit GHZ state, and a single 
M(N/M ~ [N/M J)-qubit GHZ state, giving a violation 
of 

Tr(W A p)=2 N - 1 -t N / M i. 

Note that for large M, several values give the same 
threshold. For example, all values N/2 < M < N give a 
violation of 2 N ~ 2 . It is also worth observing that since 
Wb provides a lower bound to Wa, Wb can be assigned 
the same interpretation for violating the strata of thresh- 
olds. 

Finally, one can also take a more specialized approach 
to using the witness, developing specific strategies to re- 
solve different types of entanglement. For example, were 
one to be presented with an iV-qubit pure state \ip), and 
promised that it is an Af-qubit W-state, with all other 
qubits separable, then wc might like to determine the 
value of M. This can be achieved by measuring the value 
of Wa on N — 2 different partitions. By measuring Wa on 
a subset of qubits, if that subset entirely encompasses the 
W-state, we get value Tt(WaP) — 2. However, if it only 
encompasses R of the M qubits, the value is 1 + R/M. 
Thus a systematic search using subsets of qubits 1 to n 
for 2 < n < N, the changes in value can be detected, and 
M determined. 



E. Summary 

In this section, wc have presented an entanglement wit- 
ness, Wa- Depending on the degree of violation, it de- 
tects not only full separability, but can give an upper 
bound on the number of separable qubits. The extent of 
the violation also serves to witness the presence of gen- 
uine multipartite entanglement of differing types, tuned 
most specifically to be sensitive to GHZ-likc entangle- 
ment. Since the optimization involved in calculating Wa 
is typically difficult, we presented a sub-optimal witness 
Wb that also possesses these properties, and is easier 
to calculate. Furthermore, the witness provides an up- 
per bound to the violation of WWZB inequalities - if 
Tt(Wbp) < 1 is not violated, then no WWZB inequality 
can be violated. 

In the following section, we will apply the witness to 
some simple examples. This will enable us to demon- 
strate some of the properties of the entanglement witness. 
For example, it is capable of witnessing bound entangle- 
ment. We will also be able to compare the two witnesses 
and see how tight a lower bound is provided by Wb- A 
summary of the results is depicted in Fig. [l] 
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Inverse Temperature (p) 



FIG. 1: A plot of the violation of our entanglement witness 
as a function of inverse temperature, /3, for the 5-qubit GHZ 
state (dashed) and 7-qubit cluster state (solid). Indicated in 
black are the witnesses as to the minimum number of qubits, 
{2,3,4,5}, that are entangled for values {1,2,4,8} respec- 
tively. The values of {1, 4, 8} also witness, for the GHZ state, 
the existence of {2, 3, 5}-partite entanglement. The vertical 
grey line indicates the temperature below which which either 
state can be purified by a genuine multipartite protocol. 

III. EXAMPLES 
A. Thermal Stabilizer States 

In this section of the paper we analyze entanglement in 
thermal mixtures of the stabilizer Hamiltonians with the 
help of the entanglement witnesses Wa and Wb ■ The sta- 
bilizer states, which arc the cigcnstatcs of an associated 
Hamiltonian if = —\^, n= \K n , where [K n ,K m ] = 
and Tr(K n ) = 0, enable a particularly simple descrip- 
tion of the thermal state. They are of particular inter- 
est because special cases of the stabilizer states include 
many of the important states in quantum information 
such as GHZ states, cluster states and error correcting 
codes. The examples that we present serve to illustrate 
the variety of properties that can be detected by our en- 
tanglement witness. 

The thermal state of if can be expanded due to the 
commutation of the operators, and the fact that K„ — 1, 
such that 

P = T^e-PH) = + tanh(/3/2)X»). (9) 

^ ' n— 1 

where, as usual, = fc^T, ks is the Boltzmann con- 
stant and T the temperature of the system [23] . To eval- 
uate the witnesses, we need to find the values of the Ap- 
point correlations functions and optimize over all possible 
local bases, and to achieve this we simply need to con- 
sider the products of the operators K n . 

While we have been considering calculating Wb as a 
lower bound for Wa because, in general, Wa is more 
difficult to calculate, in the following examples, it turns 



out to be no harder than the calculation of Wb- Let's 
assume that we have a minimal sequence of products of 
stabilizers that gives an A-body correlator, and perform 
local rotations such that this is the correlator of all a x s. 
Additionally, assume that all other A^-body correlators 
will come from multiplying additional terms with this 
one. The local rotations also ensure that these terms 
that we multiply by are er z s, such that the outcome is 

icr y . Thus, the overall phase is given by = cos ^| |Z|^ 

where I denotes the positions of the a z s. This cos ^| \l\j 
in Tf multiplies the identical term in the expression for 
Ag, and gives +1 since \l\ is even. Hence, all the terms are 
positive quantities, except for the term (— l) l ' k , which can 
be set positive for all I by choosing k = 0. Consequently, 
for all these cases, we can achieve 

Tr(Wp) = l T fcl 

k 

which is also an upper bound to the value of Wa, and is 
hence the optimal choice. This presents the opportunity 
of assessing how well Wb performs as a lower bound to 
Wa- In fact, given the assumption that all correlators 
are the result of the products of stabilizers, then all Tj- 
are powers of tanh(/3/2), and hence if Wa determines a 
critical temperature of Pa, then 

tanh(/3 A /2) = tanh 2 (/3 B /2). 

At the extreme of large (3a, Pb — Pa l n (2), and for small 
Pa, Pb = V^Pa, so the bound does not seem wholly 
unreasonable. 



B. GHZ Hamiltonian 

By choosing the following stabilizers, 

N 

Ki = II a W 

n=l 

K n = a^a^\ (10) 

the ground state of if is the A-qubit GHZ state (10)^ + 
)/v2. The excited states constitute the ground 
state with local operators or rji™^ 1 - 1 applied. Since 
our entanglement witness is formed from projectors on 
GHZ states, this is a natural test candidate. 

It is evident that the optimal choice of basis consists 
of the eigenstates of the Pauli operators a x and a y . This 
is particularly clear in the case of odd N, where the A^- 
body correlator necessarily includes a term due to Ki 
dln^i K n creates a correlator on all other qubits, but 
not qubit 1). Once K\ is included, multiplying by any 
other stabilizer necessarily gives an A^-body term formed 
of a x and cr y terms. All such products introduce an even 
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number of a y s so that we get 

Tt(Wap) = E ( ) tanh" +1 (/3/2) 

= tanh(/3/2) (1 + tanh(/3/2)) JV_1 
Tr(VF B( o) = tanh 2 (/3/2) (l + tanh 2 (/?/2)) ^ (11) 

In the thermodynamic limit, i.e., N — > oo, we find the 
limit of Tr(W A p) = 1 at tanh(/?/2) ~ 1/a/F. Since 
this shows that f) — ► 0, it has the interpretation that 
large systems of this form are always entangled. How- 
ever, we know from [17] that distillation is impossible if 
(3 < In (\/2 + l), which causes us to conclude that this 
persistent entanglement is bound entanglement i.e. en- 
tanglement that cannot be distilled by a multipartite 
distillation protocol. This proof arises from considering 
a particular bipartition of the qubits, and showing that 
distillation across that bipartition is impossible. How- 
ever, there exist other bipartitions across which distil- 
lation is possible. Witnesses, including Bell tests, for 
bound entanglement have previously been demonstrated 
|20l |2"T1 |2"2"] . What is perhaps most remarkable about 
this system is that entanglement becomes more persis- 
tent in larger systems. Unfortunately, the A^-body term 
K\ in the Hamiltonian is not particularly physical, so we 
should not necessarily expect to see the consequences in 
real- world systems involving local interactions. 



position. For r = 0, it is given by er 2 and a y for the end 
terms, cr x ,a y at positions 3s + 2 and cr z ,a x otherwise. 
Thus, we find that 

^|T g |=tanh r+ro (/3/2)x 

x(l + tanh 2 (/3/2)) m+r - 1 (l +tanh(/3/2)) 2 - r (12) 

and that Xf is given by the same expression, but re- 
placing all the tanh(/3/2) with tanh 2 (/3/2). In the limit 
of large TV, we find that the critical values of f3, corre- 
sponding to W A and W B read = 2.35, ^ = 1.67. 

Therefore, above the temperature Tf r — (fcs/3^.) -1 , 
the thermal cluster state does not violate two setting 
WWZB inequalities yet is entangled until at least T^. = 
(/cb/3^.) -1 , above which, the entanglement witness Wjy 
fails to detect entanglement. However, we know from 
Ref. P2] that purifiablc entanglement exists below the 
temperature given by (3 = ln(v2 + 1). Thus, for all N, 
our bound finds a lower temperature than the critical 
temperature for purification, and hence, does not detect 
any bound entanglement, in contrast to the GHZ state. 

It was experimentally shown based on the results in the 
Ref. [2] that cluster states violate some Bell inequalities 
[19] . Interestingly, these Bell inequalities are not violated 
by the GHZ state, which explains why there is a range 
of temperatures for which the thermal cluster state does 
not violate two-setting WWZB inequalities. 



C. Cluster State 

Another example of a stabilizer state is the cluster 
state. We shall restrict to the one-dimensional version, 
where the stabilizers are defined as 

K X = aWai 2) 
K N = ^ 

K _ (n-1) («) (n+1) 

From the results on purification of these states [TS] , 2D 
and 3D cluster states have exactly the same persistence 
of purifiablc entanglement, but in the present case, be- 
ing certain of having the optimal basis U is much harder. 
To compute the expectation of the witness, we consider 
three cases enumerated by r (r = Q, 1, 2). Each case cor- 
responds to the different length of the chain N r = 3m+2r 
(to is an integer). In each of these cases, the minimal 
product of stabilizers to give an A^-point correlation func- 
tion is (to + r) [25 . Further A^-body correlators can be 
constructed by multiplying by pairs of operators. For ex- 
ample, for r = 0, the basic product is K 2 K 5 K S . . . Kn-i, 
and further products can be constructed by multiplying 
terms K^ s K^ s+ \, or using the end terms K\ and K^. 
Thus the value of TV is simply a power of tanh(^), where 
the power is the number of products that have been used. 
As our Hamiltonian is defined, the best basis varies with 



D. Thermal Bose-Hubbard with Single Excitation 

Moving away from stabilizer Hamiltonians, consider a 
1-dimensional regular array of N lattice sites with pe- 
riodic boundary conditions in which we place a single 
particle (the extension to d-dimensional lattices is trivial 
because we are only using a single excitation) . This par- 
ticle is free to hop between nearest-neighbour sites with 
a constant hopping amplitude. We can define a basis |n), 
denoting that the particle is on the n th qubit. 

The cigenstates of such a system are readily expressed 

as 

\^m) = ^ W J2 e2mmn/N \ n ^ 
V n=l 

and the eigenvalues are E m — 2 cos ( ) . 

The thermal state of the system in this "position" rep- 
resentation therefore reads 

p= — ^ w) ' ( 3) 

where Z{(3) — N J2 m e~ l3Em is the partition function. 
At this stage, on each site, we associate the pres- 
ence/absence of a particle with the qubit levels |1)/|0). 
The interpretation of whether any observed violation of 
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the entanglement witness is really entanglement has been 
discussed elsewhere 23J. Instead of engaging in a com- 
plete analysis of the correlators, we can trivially observe 
that T zz „ -Z = — 1, and therefore Wb > 1- This arises 
because we know that Tr(p) = 1, and Z behaves ex- 
actly like i, except that a negative sign is introduced 
in the presence of an odd number of excitations i.e. to 
all terms. Thus, entanglement persists at all temper- 
atures. However, in contrast to the GHZ state, this 
entanglement is always purifiable, as can be proven by 
demonstrating an explicit purification protocol. Let us 
perform ^-measurements on all sites but a particular 
nearest-neighbour pair, r and r + 1, and post select on 
all measurement results being |0). This leaves the density 
matrix 



-(3E„ 



Pr,r+l — 



-2irim/N J 



2Y 



-PE„ 



The fidelity with the singlet state 



2 is 



F= i - 
r 2 



2E, 



-0E„ 



Note that the terms with negative energy have largest 
weight, and therefore F > g, and so it can be purified 
to the perfect singlet. By symmetry this can happen 
between all nearest neighbour sites, which is sufficient to 
reproduce any desired output state. 



IV. SUMMARY 

We have investigated a family Wjj of JV-qubit entan- 
glement witnesses which are functions of TV-point cor- 
relators. We have derived a lower bound Tr{Wsp) for 
Tr(Wfjp) (p is an arbitrary iV-qubit state) and demon- 
strated that the condition Tr{Wsp) = 1 separates states 
that do not violate two-setting WWZB-type Bell inequal- 
ities from the entangled states violating multi-setting 
WWZB type Bell inequalities. Both Wb and the stronger 
witness W A detect that if Tr(Wp) > 2 M+1 , at least M+l 
qubits are entangled. Similarly, above certain thresholds, 



the presence of genuine multipartite entanglement can be 
detected. 

The family of witnesses Wjj has been tested on the 
thermal stabilizer states (GHZ Hamiltonian and clus- 
ter state Hamiltonian) as well as on the one excitation 
Bose-Hubbard thermal state. In the case of the thermal 
GHZ stabilizer state, we have found that in the thermo- 
dynamic limit one can always detect entanglement with 
the help of Wjy at any finite temperature, although it is 
known that there is a finite regime of temperatures for 
which the state is entangled and purifiable. Hence, the 
state is bound entangled. Contrastingly, when we con- 
sidered the thermal cluster state, we found a critical tem- 
perature T^. below which entanglement is detected. The 
critical temperature due to the witness Wb also shows 
that if we're below T® r < t£, we have entanglement. 
However, above Tf r , no two-setting WWZB inequality 
is ever violated. In the regime between these two tem- 
peratures, the state violates multi-setting WWZB Bell 
inequalities. Finally, we have also examined the thermal 
one excitation Bose-Hubbard model, which always vio- 
lates multi-setting WWZB Bell inequalities for any finite 
temperature, although we have argued that in this case, 
the entanglement is always purifiable. 

It would be interesting to apply the entanglement wit- 
nesses Wjy to thermal states of spin Hamiltonians of 
ferromagnetic and anti-ferromagnetic type as well as to 
other models used in condensed matter physics. Another 
interesting question is whether the witnesses derived in 
the paper can be expressed as a function of the partition 
function and its higher order derivatives. Higher order 
derivatives of the partition function have not yet found 
applications in condensed matter physics. 
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APPENDIX A: DERIVATION OF W B 

In this appendix, we must prove that for our choice 
of bj: in Eqn. (|6j), we get the relation in Eqn. (|7j). It is 
convenient to write W in the Pauli basis, i.e., the basis 
consisting of elements (S> ■ • • ® &i N with tro = 1, (Ji = 
&x i &i = &y j 03 = er z . The coefficients of the expansion 
Ph...i N m this basis read 

Vh...l N =Tr((T h ®---®ai N W) = 
J^fegRe^K®.--®^^)). (Al) 

k 

Note that Re (<fe'| oi 1 <g> • • • <S> <Ji N \k)j vanishes whenever 
at least one of the indices lj = 0,3. Thus, let us define 
the vector I = (1% . . . In) such that lj = 0(1) corresponds 
to a x {(jy) at the j'th position. Thus we find the relation 



a0) = (-lf r (i)M\k') 



(A2) 



and the entanglement witnesses for a particular choice of 



(n = l,...,N), where O n is an orthogonal representa- 
tion of U n . We see that the only relevant Appoint corre- 
lation functions are those for which \l\, i.e., the number 
of a y in a-j- before the local unitary operation U, is an 
even number. 

To derive W B we put 6 £ = 2 N A £ (£fA r )~\ which still 
satisfies Eq. and is hence a sub-optimal choice of 
witness. The result is 



max (Tr(W») = max 



(A5) 

By definition Ag = 2 N - lr Ti {U^pU{Qt - Q~ )J thus 

E \h\ < ^ E Tr (u f pU(Qi + Q- k ))=2 N , (A6) 



where we have used the completeness of the operators 
, except that we remember that each is counted twice; 

k 

once for k and once for k' . Moreover, 

E4= E(-d 



If cos («)T r T„ 



l.rh.k 



2^E cos2 (f ) T r- 2N E T h ( A7 ) 



where we have made use of the formula £g(— l) fc ' ( - /+m * 1 — 
2 N d^ i7i . The notation I 6 even means summation only 

over vectors / with even Hamming weight. Combining 
the last two equations we finally get 



max 

u 



E *? . 



(A8) 



\l£even 



as required. 



U, Wjy can be written as 



^ = ^rE(- 1 ) /! "% cos (|i^^ t - 



(A3) 



Consequently, the trace of with an arbitrary density 
operator p reads 



(A4) 



The number T f = Trfjjafi^ p) is the average value of spin 
measurements along the directions given by O n x, O n y 



